Magnetic blackbody shift of hyperfine transitions for atomic clocks 
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We derive an expression for tlie magnetic blackbody shift of hyperfine transitions such as the 
cesium primary reference transition which defines the second. The shift is found to be a complicated 
function of temperature, and has a dependence only in the high-temperature limit. We also 
calculate the shift of ground-state pi/2 hyperfine transitions which have been proposed as new 
atomic clock transitions. In this case interaction with the p3/2 fine-structure multiplet may be the 
dominant effect. 
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The frequency of the ground-state hyperfine transi- 
tions used in atomic clocks (such as the cesium primary 
standard) are known to be temperature dependent 
For this reason the SI second is defined at K, and at any 
finite temperature the blackbody shift must be taken into 
account. Temperature fluctuation of the laboratory is a 
major portion of the clock error budget f5|, therefore the 
NIST-F2 cesium fountain, currently under construction, 
will be cooled to 77 K to reduce the blackbody shift. 

Recently there was some disagreement in the literature 
over the size of the electric blackbody radiation shift in 
cesium. Early measurements and ah initio calculations 
support a value about 10% higher than later measure- 
ments and semiempirical calculations (see Q for refer- 
ences). On the theory side, this seems to have been re- 
solved d, 0, in favour of the larger values. As the 
temperature of the experiment is reduced in the future 
the magnetic blackbody shift (^ T^) will become more 
important relative to the electric shift (~ T^). Hence 
this reassessment of the magnetic blackbody shift. 

In this paper we present a derivation of the mag- 
netic blackbody shift of ground-state hyperfine transi- 
tions that is valid at all temperatures (not just in the 
high-temperature limit). We calculate the effect for si/^ 
hyperfine transitions such as the 6si/2 (-F = 3 ^ 4) ^^'^Cs 
transition which defines the second (there are many other 
such clocks, including ^''Rb, ^^^Yb+, and ^^^Hg+). We 
find that the simple scaling law of the blackbody shift 
l^ujhfs ^ is only valid at high temperatures. Addition- 
ally we calculate the shift for pi/^ hyperfine transitions 
which have been proposed as clock references 0. We 
show that interaction with the P3/2 fine-structure multi- 
plet must be considered. 

The magnitude of the magnetic blackbody field is 
(atomic units Ti — e — — 1) 
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An oscillating magnetic field -B(tj) cofi{Lot) affects an 
atomic energy level via the time dependent perturbation 



V{ijj, t) = ■ B{ijj) cos{ijjt) 
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where /i is the magnetic dipole moment of the system. 
The energy is affected in the second order of perturbation 



theory (see, e.g. Q) 
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where = = fi ■ B{uj). For an atom with a sin- 
gle electron above closed shells fi = —^b{L + gsS) with 
gs = 2 and fiB = a/2 in atomic units. A general expres- 
sion for this case is presented in the Appendix. 

We first examine the case of a single si/^ orbital split 
by the hyperfine interaction with a nuclear spin /. In this 
case there are only two levels of interest, with F — /-H 1/2 
and F = / — 1/2; the next level will be separated by 
several orders-of-magnitude more than the splitting (in 
the ^^•^Cs case by a factor of 10^). Then fi = —fiB9sS, 
and one obtains 
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1^1/2 — £//_i/2 is the hyperfine splitting 
of the s-state. The total blackbody shift is obtained by 
integrating this shift over the blackbody spectrum ([T]): 
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where a = ujhfs/kT. Note that the fractional blackbody 
shift does not depend on the nuclear spin /. 

At high temperatures, a 1, the integral in ([7]) is 
analytic and one obtains 
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in agreement with Ref. [l|. By contrast, in the low tem- 
perature limit the blackbody shift has a (fcT)** depen- 
dence (and is of opposite sign): 
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FIG. 1: The integral of Eq. (|8]), normalised to the high- 
temperature limit: -^x{'^) ~ ^ J (ex_i)(ff .^^^i) • Note that 
the abscissa is 1/a: temperature increases to the right. 
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FIG. 2: Solid line: x(a) from Eq. dashed line x(a)+ax'(«) 
(used in Eq. (|9])). The high-temperature limit a ^ of both 
graphs is 7r^/6. 
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Of course, in general one can simply calculate the shift 
numerically; this can be done more easily by subtracting 
the pole. Defining F{x) — x^/{e^ — 1) the integral can 
be rewritten 



X(a) 
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since the contribution of the second term F{a) / (x^ — o?) 
is zero. x{o) is continuous and differentiable everywhere 
on the positive real axis. A graph of the integral is pre- 
sented in Fig. [1] Using Eq. ([7]) we obtain, for ^■^^Cs at 
300 K, Aujhfs/uJhfs = -1.304 x 10" 
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We now turn our attention to other clocks that use 
the hyperfine splitting of a ground pi/2-wave state as the 
reference frequency, such as those proposed in . In this 
case the blackbody radiation will again cause attraction 
(or repulsion) between the two pi/2 levels, in a similar 
fashion to the Si/2-wave case. Using the results of the 
Appendix, one obtains the first term in Eq. 

However in the py^ case there will be an additional 
shift due to interaction with the P3/2 level. In fact, there 
is the possibility of some "enhancement" of the black- 
body effect here as can be seen from Eq. ([3]): in the case 
where to ^ Ea — En = cufs the shift will be of order 
AEa ~ w/s/w^, which can be larger than the shift due to 
mixing of the pi/2 hyperfine states by a factor ujfs/ujhjs- 
However one finds that the shift for both the F = / -I- 1/2 
and F = 7 — 1/2 levels due to the P3/2 levels is identi- 
cal in second order. To go beyond second order we have 
included the differences in the energy denominators be- 
tween different hyperfine components. This affects the 
blackbody shifts at the level LUhfs/^fs, which cancels the 
enhancement mentioned earlier. One obtains for the in- 
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where ahfs = ^hfs/kT and afs — LOfs/kT. The second 
term, proportional to xio-fs) + a fsX' {'^ fs) 1 shows the ef- 
fect of the P3/2 levels on Whfs- The first part of the 
second term is due to the energy difference of the pi/2 
levels in the energy denominator, while the second part 
(~ A/iVfifs) is due to to splitting of the P3/2 levels. Here 
A is the magnetic-dipole hyperfine constant of the P3/2 
levels. Note that the electric-quadrupole terms B cancel. 

The function x'(o) in (El) arises from the expansion of 
the energy denominators and is defined by 



X{afs + ahfs) - x(a/s) 
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in the limit ahfs/afs = ujhfs/ujf. 
of x(a)+ax'(a) in Fig. H 

In the high temperature limit kT Wfs, x(a —> 

7r^/6 and x'{a ^ 0) = so from Eq. ^ we obtain 



0. We present a graph 
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Equations (|9llip show that the P3/2 levels must be in- 
cluded when calculating the magnetic blackbody shift un- 
less kT <C w/s! depending on the system considered (i.e. 
the value of 7 and A/cuhfs) it may be the dominant effect. 
In the case of Al, A/iVhfs = 0.063 [1, Q and 7 = 5/2, 
therefore the last term in pT|) is approximately 0.19. 
However in Al, cufs = 112 cm^^ = 162 K, therefore at 
300 K, afs = 0.54 and x(a) +ax'{a) ~ 0.96 « 0.59 7rV6, 
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so clearly the high-temperature limit is not appropriate. we neglect the anomalous magnetic moment (i.e. g. 
At 300 K one obtains Auj^fs/^^hfa = -2.32 x 10-^^. Then 

Our treatment of the interaction with the pa/a levels 
takes into account only the largest terms in third-order 
of perturbation theory (second-order in V and first-order 
in the hypcrfinc interaction); numerical calculation of off- 
diagonal hyperfine interaction constants is beyond the 

scope of this work. In any case usually the non-diagonal AEa — 7 \ {bmb \ fJ-B{J + S) ■ B\ 

hyperfine matrix element {pi/2 \Hhfs\p3/2) is significantly 
smaller than the diagonal one {tUhfs for pi/2), therefore we 
do not expect the result to change significantly. However 
when a clock is produced a more accurate third-order 
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calculation will be necessary. In the meantime the last _ Mb-^ i'^) \^ — 

term of ([9]) (~ A/ujhfs) may be considered an estimate of 6 ^ ° {Ea — E^Y ~ 

the error. 
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APPENDIX: MAGNETIC BLACKBODY SHIFT 

The total magnetic blackbody shift of a level \ama) is 
given by equation ([3]). For electrons /x = —^ib{J + S) if 



Consider the interaction of levels in an atom with nuclear 
spin /. If we denote the angular quantum numbers of \a) 
with L, J and F, and those of \h) with L', J' and F\ 
then 



2) 



where P = J + L + + F + F' + 21 + I and the square brackets [J] = (2 J + 1) 
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